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Abstract 

e. 

^ ■ Using density functional theory, we investigate fluctuations of the ground-state en- 

ergy of spin-polarized, disordered quantum dots in the metallic regime. To compare 

O ■ 

to experiment, we evaluate the distribution of addition energies and find a convolution 
of the Wigner-Dyson distribution, expected for noninteracting electrons, with a nar- 
rower Gaussian distribution due to interactions. The third moment of the total dis- 
tribution is independent of interactions, and so is predicted to decrease by a factor of 
(2 — 57r/8) / (2 — 6/71") ~ 0.405 upon application of a magnetic field, which transforms from 
the Gaussian orthogonal to the Gaussian unitary ensemble. 



PACS numbers: 73.61.-r, 73.23.Hk, 71.15.Mb 



1 



The interplay of disorder and electron-electron interactions in quantum dots has re- 
cently attracted much attention. Experiments using quantum dots - small islands fabri- 
cated in a two-dimensional electron gas [IJ - measure the spacings between conductance 
peaks in the Coulomb blockade region. Since the peak spacings reflect differences between 
ground-state energies for different numbers of electrons one cannot apply random matrix 
theory |2[ to evaluate the spectrum of peak-spacing fluctuations. Indeed, experiments 
find a more symmetric distribution than the Wigner-Dyson form. Experiments disagree, 
however, on the magnitude of the fluctuations || |, EL H]. Sivan et al. ||, observed 
fluctuations several times as large as the inferred mean level spacing ( A ) , and concluded 
that the fluctuations are a fixed percentage 10 — 15% of the total charging energy e 2 /C, 
where C is the dot capacitance. Similar results have recently been obtained by Simmel 
et al. [H in small Si dots. In contrast, Patel et al. [f| found fluctuations in GaAs dots 
comparable to the mean level spacing. 

Theoretical treatments also disagree regarding the magnitude of the peak-spacing 
fluctuations. Sivan et al. [|3| found large fluctuations scaling as (0.10 — 0.17)e 2 /C for 
a small lattice model. Similar results were found by Koulakov et al. for the classical, 
strong interaction regime r s ^> 1 @ where electrons form a Wigner lattice f| . Blanter et 
al. used the random phase approximation (RPA) |10| for weakly interacting dots and 
concluded that, for dimensionless conductance g 3> 1 flllj], the contribution to fluctuations 
from interactions should be parametrically smaller than the mean level spacing (A ). 
While the above results can be reconciled as applying to different regimes of r s and g, 



recent work employing the self-consistent Hartree-Fock equations [JT2], [13], |Tj] found peak- 
spacing fluctuations several times as large as (A ) even for r s ~ 1 and g ^> 1 where RPA 
should still provide a good approximation JTO) . 

The purpose of the present article is to clarify the origin, magnitude, and distribu- 
tion of peak-spacing fluctuations in spin-polarized disordered quantum dots in the regime 
g ~ 1 and r s ~ 1. Density functional theory (DFT) provides us with accurate ground- 
state energies including electron-electron interaction, confinement, and disorder for real- 
istic quantum dots. We find that the distribution of peak spacings is the convolution 
of a Wigner-Dyson distribution, expected for noninteracting electrons, with a narrower 
Gaussian distribution due to interactions. The width of the Gaussian is accurately given 
by the fluctuations in the screened Coulomb interaction between a pair of electrons at 
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the Fermi energy - a result which will also apply to unpolarized quantum dots. The 
total peak spacing fluctuations are hence smaller than the mean level spacing (A ). Use 
of an unscreened interaction between electrons, either direct or exchange, is found to 
greatly overestimate the magnitude of the fluctuations. Furthermore, since interactions 
add a symmetric contribution to the distribution of peak-spacing fluctuations, the third 
moment of the total distribution is independent of interactions. Hence, we predict that 
experimental application of a magnetic field will reduce the third moment by a universal 
factor of 0.405, corresponding to a change from the Gaussian orthogonal ensemble to the 
Gaussian unitary ensemble. 

The ground-state energies of spin-polarized, disordered quantum dots are obtained 
within density functional theory with the exchange-correlation part of the electron-electron 
interactions treated in the local-density approximation. Specifically, we solve the follow- 
ing Kohn-Sham equations [l5[ numerically, and iterate until self-consistent solutions are 
obtained 111 



» 2 v ^/^ r < + ^W; 



2m* k J |r — r'| <5p( r ) 

where the density is 



^(r) = e i * i (r), (1) 



N 

p(r)=Y,\Mr)\ 2 - (2) 



Here E xc [p, (] is the exchange-correlation energy functional JT7| with local spin polarization 



C(r) = 1. The summation in the density (0) is taken over the N lowest energy Kohn- 
Sham orbitals. In previous work |Tj|, we have shown that the DFT method gives very 
accurate ground state energies for clean parabolic GaAs quantum dots, in agreement with 
exact calculations for up to five electrons ||T^|. Comparison with quantum Monte Carlo 
calculations [2U| confirms that DFT is valid for interaction strengths up to (e 2 / ntoj/TwoQ = 
6 (r s ~ 8) and up to N = 8 electrons. 

The external potential for our disordered dots is the sum of a confining parabola and 
multiple "impurity" potentials each with a Gaussian profile: 

1 1 JVimp / |r-r| 2 \ 
Kxt(r) = -m^y + — £ 7l ■ exp ^ _ j . (3) 

The impurity potentials are randomly distributed with density n- irav> = 1.03 x 10~ 3 nm -2 
and strength 7$ uniformly distributed on [— W/2, W/2] with W = 10h 2 /m*. The width 



of each impurity is taken as A = £q/{2\/2), where £o = yh/mFujQ ~ 19.5nm. Here we 
use the effective mass for GaAs, m* = 0.067m, and hujQ = 3.0meV. The strength of the 
Coulomb interaction is controlled by changing the dielectric constant k, where k = 12.9 
for GaAs. The resulting dimensionless interaction strength is measured by (e 2 /n£o)/tuJo 
or r s (= 1/ \jirpoOg) where a* B = ti 2 n/m*e 2 is the effective Bohr radius, and p is the 
electron density at the center of the dot. From a scattering phase-shift analysis we find 
the mean free path of electrons I = VpT ~ 170nm to be slightly larger than the dot 
diameter L = 120 — 160nm, where the dot diameter increases with r s . Therefore the dots 
are marginally in the ballistic regime and have a dimensionless conductance g = 2 — 4 

[n|. 

At low temperatures, electron hopping into a dot containing AT — 1 electrons is sup- 
pressed except when the ground-state free energy ^(A^ — 1) — (N — is equal to the 
ground state free energy for N electrons E(N) — N/j,. This degeneracy condition deter- 
mines the position of the A^th conductance peak as a function of the electron chemical 
potential hn = E(N) — E(N — 1), or equivalently, as a function of an applied gate voltage 
|[21|| . The increase in \x needed to put an extra electron in the dot, which we will refer to 
as the addition energy A, is given by A = E(N + 1) — 2E(N) + E(N — 1). From our 
solution of the Kohn-Sham equations, the ground-state energy of a dot with A^ electrons 
is obtained from 

W-E,-£/f^-/*)^U + *- (4) 

We consider fluctuations of the addition energy for A^ = 10 electrons. Thus for each 
realization of disorder we calculate A = E(ll) — 2E(10) + E(9). The disorder average is 
taken over more than 1,000 different impurity configurations. As a check of accuracy, we 
have confirmed that the ground-state energies obtained from DFT for disordered quantum 
dots with N = 2 and 3 are in good agreement with exact diagonalization results for 
< r s < 5. 

Figure 1 shows the distribution of addition energies A for interacting dots with 
(e 2 / 'k£o) /hujQ = 2.39 (right), and the distribution of A for noninteracting dots of the 
same size (left). In the inset, we show the charge density p(r) for one realization of 
disorder. The distribution of level spacings Ao in the noninteracting dots has the Wigner- 
Dyson form, while that in the interacting dots is somewhat more symmetrical. The 
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symmetry continues to increase with increasing interaction strength. While interactions 
considerably enhance the average addition energy (A) ~ 6.50(A ), the fluctuation in 
the interacting case 5A = \J (A 2 ) — (A) 2 is only ~ 13% larger than the noninteracting 
fluctuation 5A = ^/(A§) - (A ) 2 . 

Figure 2 shows the average addition energy (A), its rms fluctuations 5 A, and its 
third moment for disordered dots as a function of the Coulomb interaction strength 
(e 2 / k£ ) /Tuoq. For comparison, we have also plotted results for disordered, noninteracting 
dots of the same size, which we obtain as follows. First we find the effective potential 
for the clean dot, without impurities. Then we solve for the single-particle level energies 
e° for this effective potential plus the random impurity potentials in (El). The addition 
energy is simply given by Ao = ejv+i — with N = 10. All dot sizes satisfy the relation 
5Aq = ^Ja/ti — 1 (Ao) — 0.52(A ) predicted by random matrix theory for the Gaussian 
orthogonal ensemble 0. In Fig. 2(a), the average addition energy in the noninteract- 
ing case is seen to decrease with increasing Coulomb interaction. This is because the 
increasing Coulomb repulsion among electrons causes the dot to grow and hence the level 
spacing to shrink. The average addition energy in the interacting case increases consid- 
erably with Coulomb interaction strength, as expected from the classical electrostatics 
e 2 

relation (A) ~ — where C is the capacitance of the dot. However, Fig. 2(b) shows that 
the interactions only slightly increase the addition-energy fluctuations. For GaAs, r s ~ 2, 
the enhancement is only about 10%, in rough agreement with the experiment of Patel et 
al. §. 

To understand the magnitude of addition-energy fluctuations, we use the phenomeno- 
logical framework presented by Blanter, Mirlin, and Muzykantskii || for the regime 
r s < 1, where RPA is valid, and show that it applies to the DFT results at least up 
to r s ~ 5. Consider first a dot containing N — 1 electrons in the ground state. Addition 
of the iVth electron to form the N electron ground state requires an electron chemical 
potential fi^. To form instead the first excited state of N electrons requires the higher 
chemical potential /i^ + Ae. For an ensemble of disordered metallic dots, Ae will have 
Wigner-Dyson statistics, with (Ae) equal to the mean noninteracting level spacing (A ), 
since the lowest excitation of a Fermi liquid is a single electron promoted across the Fermi 
surface. 

The addition energy A is the increase in chemical potential from fi N required to add 
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one more electron to the dot and thus form the N+ 1 electron ground state. This (iV+ l)st 
electron must have an extra energy Ae to occupy the lowest empty level plus an extra 
energy Un,n+i due to its Coulomb interaction with the Nth electron. The total addition 
energy A will be approximately given by the sum of these two contributions, 

A~ Ae + U N , N+1 . (5) 

The distribution of Ae is given by the Wigner-Dyson distribution of level spacings for a 
noninteracting dot of the same size. The average interaction energy (Un,n+i) is the ca- 
pacitive charging energy e 2 /C. We estimate the fluctuations in Un,n+i by calculating the 
screened Coulomb interaction between two electrons at the Fermi surface 0. Specifically, 
we treat the screening effect in the Thomas-Fermi approximation as 

Ux F N+x = ej ip N {r)p% +l {v)dr. (6) 



The screened potential due to the Nth electron in Fourier representation is [22 



2?re p° N (q) 
K |q| + g 

where p%(r) = |0^(r)| 2 is the density of the iVth single-particle wavefunction 4>%(r) of a 
noninteracting disordered dot. The Thomas- Fermi wavevector is go — (27re 2 /K) (dn/dp) = 



l/a* B . It is found that the fluctuation 5U TF = J ((^vjv+i) 2 ) — (^vjv+i) 2 is always con- 
siderably smaller than the noninteracting level-spacing fluctuation 5A up to at least 
r s ~ 5. The total fluctuation estimated as 5A TF = J (<5A ) 2 + (5U TF ) 2 is shown in Fig. 
2(b) by crosses. We see that the fluctuations in the Thomas- Fermi screening model agree 
well with the DFT results with no free parameters. This supports the picture || that 
the addition-energy fluctuation arises from two quasi-particles above a filled Fermi sea 
interacting via a screened Coulomb potential. 

Within this picture, the increase of the fluctuation of [/jv,jv+i with increasing interac- 
tion strength leads naturally to greater symmetry of the distribution of addition energies. 
Numerically, we find that that the distribution of i/jviv+i nas a symmetric Gaussian 
form. Hence, in agreement with Eq. (|J), we observe that the addition-energy distribution 
function .P(A) is always extremely well described by the convolution of a Wigner-Dyson 
distribution for level spacings Ae, 

ii Ae 



with a Gaussian distribution for interaction energies Unn+i, 



,2 



I (C/jV,iV+l-(A) + <A ))^ 

P Ga UU N ,N+i) = -7=7- e . (9) 

V 2,no U 



The result for the distribution of addition energies is 



A) 



1 /7T <?£/ -S A 2 .. f A2 

.g 4 Q (A >2 J g 2a(«E7)2 _)_ / 



1 + erf( 



2a5CT 



(10) 



2V2a(A ) 2 [ V2a5[7 

Here a = ttSU 2 / (2(A ) 2 ) + 1 and A = A - (A) + (A ), where (A) is the center of the 
distribution and SU is a fitting parameter giving the width of the fluctuations of £/jv,jv+i- in 
the noninteracting case, SU = so that -P(A) = Pwd(A) as expected. In the other limit, 
P (A) becomes nearly symmetric for sufficiently large SU. In Fig. 1, we show P(A) given 
by Eq. (pTDj ) as a dashed line. It is seen that the DFT distribution is described very well by 
([10|) with the best fit value of SU = 0.13meV very close to the value SU^n+i = O.lOmeV 
estimated from the Thomas-Fermi screened Coulomb interaction between two electrons 
at the Fermi surface. 

To test whether the distribution of addition energies is well described by the sum of 

noninteracting level spacings and a symmetric distribution due to interactions, we propose 

to compare the third moment of the distribution P(A) with and without a magnetic 

field B± normal to the plane of the dot. Since the interaction part, coming from the 

screened Coulomb interaction in our picture, is symmetric it does not contribute to the 

third moment of P(A). Therefore, the ratio -77— -. \ \1}, B± ^° should take the value 

((A — (A)) 3 ) B±=Q 

(2 — 57r/8)/(2 — 6/7r) ~ 0.405 which applies to level spacings taken from a Gaussian 
orthogonal ensemble (B± = 0) and a Gaussian unitary ensemble (B± ^ 0) 0. Since our 
results apply only to the case of spin-polarized electrons, it is necessary to apply a large 



magnetic field in the plane of the dots, or to spin polarize the nuclei [23|. The result can 
also be tested numerically, e.g. by exact diagonalization studies as in Ref. 0. 

Existing diagonalization studies for spin polarized electrons on small lattices find 
addition-energy fluctuations SA ~ 0.15e 2 /C ||. For comparison, the Coulomb contri- 
bution to the fluctuations found by DFT are much smaller SA ~ 0.03e 2 /C at r s ~ 2. 
This difference may be attributed to differences in the strength of disorder: while the 
dimensionless conductance in our dots is g = 2 — 4, we estimate g = 0.1 — 0.3 in Ref. 0. 
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Various theoretical estimates give fluctuations SA oc (Ao)/g[pl| or 5 A oc (A ) / \fg ||. In 
either case, the discrepancy between exact diagonalization and DFT can be attributed 
to the order of magnitude difference in the dimensionless conductance g in the samples 
studied. The experiments on GaAs []3], f|, |5| have r s ~ 1 and g > 1, and thus fall more 
closely in the range of interaction strengths and dimensionless conductance treated in this 
paper. 



Recently, several calculations fll|, [13| [L4[ based on the self-consistent Hartree-Fock 
(SCHF) equations have found large fluctuations, up to 5 A ~ 0.2e 2 /C, in the same range 
of r s we consider. In one case JT3|, the dimensionless conductance is estimated to be 
g ^> 1, a regime where DFT predicts fluctuations an order of magnitude smaller. As 



pointed out by Walker, Montambaux, and Gefen ||12|| , the exchange interaction in the 
SCHF equations is unscreened. To test whether the lack of exchange screening in the 
SCHF approach could be responsible for the discrepancy with DFT, we have calculated 
the unscreened exchange interaction between two electrons near the Fermi surface in our 
dots 

In Fig. 2(b), we have plotted as open circles the fluctuations taken by summing the 
unscreened exchange interaction ( [TTj) with the noninteracting level spacing [25]. It is 
clear that for r s > 1, the unscreened exchange interaction noticeably overestimates the 
addition-energy fluctuations. In contrast, density functional theory correctly accounts 
for screening within the electron gas, including exchange interactions ||26|| . These results 
suggest that the unscreened exchange interaction in the SCHF approach may generally 
lead to an overestimate of the addition-energy fluctuations. 

In this work, we have neglected external screening by gates or electrodes. This simpli- 
fication should be valid as long as the distance to external conductors is larger than the 
diameter of the dot. In the opposite limit, it is essential to consider external screening, 
but this may be done by a simple modification of the 1/r potential between electrons. 

In conclusion, we have studied the electronic states of spin-polarized, disordered quan- 
tum dots using density functional theory and investigated the fluctuation of the ground- 
state energies. We have found that electron-electron interactions increase the fluctuation 
of addition energies by no more than 25%, up to r s ~ 5, even though the average addition 
energy is increased by a factor of 10. The addition energy is well approximated as the 
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sum of the noninteracting level spacing and the screened Coulomb interaction between 
two electrons at the Fermi surface. Hence the distribution of addition energies is the 
convolution of a Wigner-Dyson distribution of level spacings with a Gaussian distribution 
of interaction energies. Since the latter is symmetric, it does not contribute to the third 
moment of the addition-energy distribution. The third moment is therefore predicted to 
decrease by a universal factor of 0.405 on application of a magnetic field which transforms 
the dot from the Gaussian orthogonal to the Gaussian unitary ensemble. For quantum 
dots having larger numbers of electrons, whether spin-polarized or not, we anticipate that 
the decrease of the screened Coulomb interaction-energy fluctuations occurs as fast as the 
decrease of the level spacing and thus the present results are also applicable. 

We acknowledge I. L. Aleiner, B. L. Altshuler, R. Berkovits, C. M. Marcus, and 
M. Stopa for useful comments and suggestions. 
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Figure Captions 

Figure 1: Distribution of addition energies A to add the 11th electron for interacting dots 
with (e 2 / kIq) I '%ooq = 2.39 (right) and A for noninteracting dots of the same size 
(left). The energy binwidth is 0.05 [meV]. The dashed lines show the distribution 
function obtained from Eq. (10). Inset: The charge density profile p(r) for N = 10 
electrons with one configuration of impurities. 

Figure 2: (a) Average addition energy (A), (b) fluctuation 8 A, and (c) cube-rooted third mo- 
ment \J ((A — (A)) 3 ), as a function of electron-electron interaction strength (e 2 / Ki )/hu . 
The measure of interaction strength r s (— l/^/np a* B ) indicated by arrows in (a) 
can be applied to the data in all panels. For each data point, the disorder average 
is taken over more than 1,000 different impurity configurations. At each r s , the 
noninteracting data are taken for dots of the same size as the interacting dots, and 
the relation 8A ~ 0.52(A ) expected for noninteracting level-spacing statistics is 
always satisfied. Also plotted in (b) are the fluctuations due to noninteracting level 
spacings plus the screened Coulomb interaction between two electrons at the Fermi 
surface (crosses), or the unscreened exchange interaction between the two electrons 
(open circles). 



13 




Addition Energy [meV] 




12 3 4 

(e 2 M>)/# 



